Bounds on Localized Modes in the Crystal Impurity Problem 

Harry J. LipkirQ 

Department of Particle Physics, Weizmann Institute of Science, Rehovot 76100, Israel 



Philip D. Mannheir 

Department of Physics, University of Connecticut, Storrs, CT 06269 

(Dated: February 7, 2006) 

Abstract 

Using general properties of the crystal site representation normal mode matrix, we provide some 
very simple bounds on localized modes in simple, body-centered and face-centered cubic crystals 
with substitutional point defects. We derive a trace condition constraint on the net change in crystal 
eigenfrequencies caused by the introduction of a defect, with the condition being a completely 
general one which holds for any combination of central and non-central crystal force-constants 
and for all- neighbor interactions. Using this condition we show that the sufficient condition for 
producing localized modes in an arbitrary cubic crystal by a mass change at the defect site is that 
the defect mass be less than one half of that of the host atom mass which it replaces, and that the 
sufficient condition for producing localized modes in an arbitrary cubic crystal by force-constant 
changes alone is that the defect site self force-constant be greater than twice that of the pure crystal 
self force-constant of the host atom which it replaces. 
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I. INTRODUCTION 



The substitutional insertion of a point defect impurity into an otherwise perfect host 
crystal will typically modify the spectrum of the 3N normal modes of the crystal,— leading 
in certain circumstances (such as the insertion of impurities which are lighter in mass than 
the host atoms which they replace or which are more strongly coupled to the host crystal 
atoms than the ones they replace) to the generation of modes with frequencies which lie 
beyond the band maximum u ma3L of the crystal. Such modes will not be plane waves which 
propagate throughout the crystal, but will instead fall off exponentially fast away from the 
defect site and thus be localized to it. Moreover, with the rest of the crystal atoms not 
participating appreciably in such localized modes, the intensity of the defect in such modes 
will be N times larger than the intensity it would otherwise have had in a crystal plane wave 
mode, to thus give the localized mode enough intensity to render it observable. While such 
modes could be of relevance for phenomena such as the Mossbauer effect associated with 
the insertion of Mossbauer active defects into host crystals, historical recoil-free fraction 
Mossbauer studies only involved an averaging over all the lattice modes of the system, to 
thereby only allow one to infer the possible presence of localized modes indirectly. However, 
with the advent of dedicated synchrotron rings it became possible to monitor Mossbauer 
active systems mode by mode directly; and via nuclear resonant inelastic x-ray scattering 
studies, modes lying beyond the host crystal phonon band maximum have now explicitly 
been seen in the Fe/Cu and the 57 Fe/NiAl systems.— 1 ^ Consequently, knowing whether the 
insertion of a defect into pure crystal hosts might generate localized modes can be of great 
value for such studies. 

In this paper we use a very straightforward trace technique to enable us to identify some 
very general conditions under which such localized modes can be produced. We shall restrict 
our study to the most straightforward case of single point defects which are substitutional^ 
inserted in the three primitive simple, body-centered and face-centered cubic crystals (these 
specific cases being amongst the most commonly experimentally studied ones), though with 
our approach being quite generic, it could in principle be adapted to encompass other crystal 
structures as well if desired. We shall provide results for various combinations of mass and 
force-constant changes, some which are specific to nearest-neighbor force-constants and some 
of which are general to all-neighbor force-constants. In Sec. II we derive the trace condition 
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on which all of our results are based, with the provided relation (Eq. (j!5|) below) being an 
exact, all-neighbor relation which permits an arbitrary mass change (M replaced by M') at 
the defect site and arbitrary force- const ant changes between the defect atom and any other 
atom in the entire crystal, despite which the resulting relation only involves the self-force 
constant A xx (0,0) at the defect site and the change ^4^(0,0) in it. In Sec. Ill we apply 
our trace condition to the case of an isotopic substitution where the only change is that 
in the mass at the defect site, to show that while having a lighter impurity (M' < M) is 
necessary for the generation of a localized mode, it is the condition M' < M/2 which is 
the sufficient one. In Sees. IV and V we apply our trace condition to some straightforward 
nearest-neighbor force-constant change cases for which there are extremely simple exact 
analytic solutions to the impurity problem (some typical examples of which being crystals 
with force-constants which are central or which are isotropic), to find that for all of them 
the sufficient condition for localized mode production by force-constant changes alone is 
given as the requirement that A' xx (0, 0)/A xx (0, 0) be greater than 3/2. Finally in Sec. V 
we examine more complicated force-constant change cases and go beyond nearest-neighbor 
force-constants, and while these situations do not admit of as straightforward a treatment 
as the cases considered in Sec. IV, nonetheless for them we are still able to extract a 
general sufficiency condition for the generation of localized modes in an arbitrary crystal by 
force-constant changes alone, namely that A' xx (0, 0)/A xx (0, 0) be greater than two. 



II. DERIVATION OF THE TRACE CONDITION 



Since the very introduction of a defect breaks the translation invariance of the lattice, 
to derive a condition such as one involving a trace, we will need to work entirely in the 
coordinate space crystal site representation. To actually derive the trace condition we recall 
that in the harmonic approximation the equations of motion for the displacements from 
equilibrium e~ ujt u a (£) of the atoms of a pure 3N-dimensional cubic crystal lattice are given 
by 

Y, [Aapit, tj - w 2 M(£')5 a(3 5(£, £')} u p {£') = , (1) 

where £ ranges from to N — 1, a — x, y, z, M(£) is the mass of the atom at site £, and 
A a p (£,£') are the second order force-constants. Since Eq. (0) is an eigenvalue equation for 
the frequencies u 2 , on defining a matrix B a p(£,£') = M(£)8 a p8(£,£ f ), we immediately see 



that the sum of the eigenfrequencies of the crystal is given by 

3iV 

= TrAB- 1 . (2) 

8=1 

For a pure crystal in which all masses are equal to a common M(£) = M and all self force- 
constants are equal to a common A xx (£,£) = A yy (£,£) = A zz (£,£) = A xx (0,0), Eq. (J2J) 
reduces to 

— v 2 = ^(o 1 o) 

3iV M K ' 

With the left-hand side of Eq. (J3J) being recognized as /X2, the second moment of the density 
of states, Eq. (J3J) thus recovers the well-known relation for pure crystals 4 

For the system with a substitutional point impurity of mass M' located at the origin of 
coordinates and changed force- const ants A' a/3 (£,£'), the displacements from equilibrium are 
now given as e~' lul ' t u a (£), with Eq. (JIJ) then being replaced by 

£ [A a/3 {£, £') - w' 2 M6 aP 6(£, £')] u p {£') = £ V a/3 {£, £')u^£') , (5) 
p,e' /3, e 1 

where the changes from the pure crystal case are described by the perturbation 

V«p(e, £') = -w' 2 (M - M')5 aP 6(£, 0)6(£', 0) + A aP (£, £') - A' af3 (£, £') . (6) 

On now defining a matrix B' ap (£, £') = M5 a p5(£, £') + (M' - M)5 a p5(£, 0)5(£', 0), the sum of 
the eigenfrequencies of the perturbed crystal is then given by 

3N 

$>? = TrA'[B')- 1 , (7) 

i=i 

with the change in the sum of the eigenfrequencies thus being given by 

3N 

YM-ufi = TrA'[B'\- 1 -TrAB- 1 . (8) 

i=l 

The utility of Eq. (jSJ) is that its left-hand side measures whether modes are shifted to 
higher or lower frequency, while on its right-hand side it is only the sites which are explicitly 
involved in V a p(£, £') which do not drop out of the difference between the two 3N-dimensional 
traces. Our task is thus to first evaluate the right-hand side of Eq. (jHJ) and then to seek 
constraints on its left-hand side. 
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In the simplest case where the only change is a mass change at the defect site, it is only 
the defect site contribution itself which is not cancelled in the trace difference, with the full 
change in frequency immediately being given by 

, 2 2 _ 3A XX {0,0) 3^(0,0) 

i=l 



3N 

£ 



(9) 



M< M 

When there is also a change in force- const ant at the defect site, we recall that by Newton's 
third law of motion there must also be changes in the force- const ants at other sites too. 
Moreover, this same Newtonian law entails that the force-constants of a perturbed system 
with a defect have to obey 

A' a p(£,£) = -J2A' a e(£',£) (10) 
for all i in exactly the same way as the pure crystal force-constants have to obey 

A aP (£,£) = -J2A a p(£',£) . (11) 

With the x, y and z directions being equivalent in cubic crystals, in the general case which 
allows for arbitrary force-constant changes, Eq. (jHj) takes the form 



3JV 



A' xx (0,0) A xx (0,0) 



i=l 



M> 



M 



A' xx (l, 1) A xx (l, 1) + 4^(2,2) _ 4^(2,2) + 



M M M M 

As regards the pure and the perturbed force- const ants, we recall that in terms of the two- 
body interatomic potential <j)(r), the force-constants between an atom vibrating around site 
R a {£) and one vibrating around the origin are defined as 



A a p(£,0) 



D 2 



du a {£)du p {£) 



0"(r) <f)'(r) 



R a (£)Rp(£) - ^-Sap , (13) 



as calculated at the equilibrium separation r. Then, since the effect of the introduction of 
the defect is to alter the two-body potential between the defect and the host atoms while 
not affecting the potentials between any two host atoms themselves, the only non-self force- 
constants which will change will then be the A' af3 (£, 0) with I ^ 0, with Eq. (jl())l then 
obliging the self force-constants at the defect site and those at the I ^ sites to respectively 
change as 



4^(0,0) 



A' 



-4*&0)- £ A aP {£,£') 



-A' a p(£,0)+A aP (£,£) + A aP (£,0) . (14) 



Given these relations and Eqs. (JTUj) and Eq. (|T2")) can thus be simplified to 

alm 2 ^(°'°) ri5 ^ 

Equation (|15|) is our main result, and is derived here with no restriction at all on the number 
of neighbors of the defect for which the force-constants might change. Nor does it presuppose 
any relation between <fr"(r) and <fr'(r)/r. Despite the fact that Eq. (fTB^) conveniently only 
involves the self force-constant at the defect site, nonetheless it is an all-neighbor result, one 
which additionally holds for any combination of central and non-central force-constants. 5 



III. APPLICATION TO THE MASS DEFECT CASE 

In the treatment of crystal impurity problem it is conventional to solve Eq. (J5J) by the 
lattice Green's function technique, and since we will use some of its aspects to constrain 
the left-hand side of Eq. (|15|). we briefly recall the procedure. One first introduces the 
dynamical matrix of the pure crystal 

^(k) = ^E^(0^)e~ iCliW (16) 

as expressed in terms of the phonon modes k of the translational invariant pure crystal, and 
then defines its eigenvectors and eigenvalues according to 

££Uk)4(k) = ^ 2 (kK(k) , 

p 

E< J (kK'(k) = 5 jf , £^'(k)^(k) = <W , (17) 

a j 

and then uses these eigenvectors and eigenvalues to construct the pure crystal lattice Green's 
functions according to 

M";^) = ^X; ^—^ ] < 18 > 

as summed over the polarizations j = (1, 2, 3) and momenta k of all the modes in the 
Brillouin zone. As constructed these Green's functions obey 



^^(O,^; = MLU 2 g a , a (uj;0,£') + J^L , (19) 



and thus immediately allow us to solve Eq. (J3J) in the form 

«aW= E gap(u , ]e,e')Vf ly (e , ,e")u J (e") , (20) 

with the eigenmodes being given as the solutions to the (3iV-dimensional) determinantal 
condition |1 — Go^l = as written in an obvious notation. 

For the simplest case of just a change in mass at the defect site, Eq. (120)) requires the 
three components of the defect displacement vector to obey 

u x (0) = -(M -M')u ,2 g xx (u';0,0)u x (0) , 
Uy(0) = -(AT - MV\y(^;0,0K(0) , 

u z (0) = -(M-M')u ,2 g zz (u'; 0,0)^(0) , (21) 

where the cubic symmetry of the host lattice requires that g xx (u'; 0, 0), g yy (ui'; 0, 0) and 
g zz (uj'; 0, 0) all be equal to each other, and thus that each one of them be given as 

g xx (u'; o, 0) = — £ ^ ^ J = - [fc(u/; 0, 0) + g yy (u'; 0, 0) + g zz (u'; 0, 0)] 
1 [gggH(k ) + <(k)^(k) + a^(k)^(k)] 



3iVM^ h 2 (k)-^' 2 ] 

— ^E — ^ 

3iVM i? h 2 ( k ) - ^ 

1 E^ = 4f^V^ (22) 



^ (w? - uj' 2 ) M Jo (to 2 - to 12 ) 

as now summed over N threefold degenerate pure crystal eigenmodes io> 2 .— And with the de 
terminantal condition then being given by |1 — GqV\ = (M' /M) 3 [l — p(to' 2 )S(uj' 2 )] 3 (l) 3N ~ 3 = 
0, in the mass defect case the perturbed simple, body-centered and face-centered cubic crys 
tal modes are thus given as the solutions to the familial 



where p(to' 2 ) is given by 
and S(ui 12 ) is given by 



1 - p{u' 2 )S{uo' 2 ) = , (23) 



P^' 2 ) = mJ~ 1 ( 24 ) 



1 N 2 

S(u' 2 ) = -1 - Mto' 2 g xx (to'; 0, 0) = - E , /2 1 2 , , (25) 



FIG. 1: Plot of S(lo' 2 ) versus uj' 2 . The line AB corresponds to the value of l/p(u)' 2 ) in a typical 
lighter defect case with M' < M, while the line CD corresponds to the l/p(u' 2 ) associated with 
an infinitely heavy defect case. 

as summed over the N eigenfrequencies u 2 of the pure crystal. 

For our purposes here we note that the function S(uj' 2 ) is divergent at every pure crystal 
frequency, to thus contain the set of N asymptotes indicated in Fig. (1). S(uj' 2 ) starts out 
equal to minus one at uj' 2 = 0, and decreases to minus infinity at the first asymptote, with 
S{ui' 2 ) thus having no zero below the lowest pure crystal band mode. Since the derivative 
of S{u' 2 ), viz. 

jgfeg) = _lf (26 ) 

is negative definite, the function S^u/ 2 ) falls monotonically between asymptotes, to thus 
have a zero somewhere on the OX axis between any two adjacent pure crystal modes, while 



8 



beyond the last pure crystal frequency eigenmode at u^ax the function S(uj' 2 ) is positive, 
falling to zero at uj' 2 = oo. As introduced, the function S(uj' 2 ) thus has N — 1 zeroes inside 
the pure crystal phonon band and one zero outside. 

In consequence of this structure for S(u' 2 ), we see that the curve 1/ p(uj' 2 ) = M'/ (M — M') 
will intercept the S(uj' 2 ) curve N times no matter what the value of M'/Mp with Eq. 
(12 1|) thus always having the same number of mode solutions as the pure crystal (viz. N 
solutions for each of the three u x (0), u y (0) and u z (0) sectors of Eq. (j2U)), with Eq. (|2"3"j) 
and the threefold degeneracy of Eq. (J2Tj) thus accounting for the complete spectrum of 
eigenfrequencies of the perturbed crystal. When M' is less than M (typical intercept AB) 
the N modes which satisfy Eq. (J23|) are all shifted to higher frequency compared to their pure 
crystal counterparts, with both sides of Eq. (jUJ) thus consistently being positive. Similarly, 
when M' is greater than M (typical intercept CD) the N modes which satisfy Eq. are 
all shifted to lower frequency compared to their pure crystal counterparts, with both sides 
of Eq. (JHJ) then consistently being negative. Finally, once the eigenfrequencies have been 
found, a return to Eq. (JHJ) will then allow a determination of all of the 3N displacements 
u a (£) in every mode, to then completely specify the crystal displacements of the perturbed 
crystal. 

Since the zeroes of S(uj' 2 ) other than the final one at uj' 2 = oo lie between adjacent 
asymptotes, an initial glance at Fig. (1) would suggest that when shifted upwards, each of 
(N — 1) modes other than the mode at the band maximum itself could be shifted upward as 
far as the immediate next asymptote, to suggest that the quantity Y.fji 1 ^? ~ uf] could be 
as large as u; 2 ^ — uj\ ~ u;^ ax (the net shift between the highest and lowest in-band frequency 
modes). However, we will now show that in fact the net shift must be much less than this, 
an outcome that will sharply constrain the implications of Eq. (j!5|) . To this end we need 
to obtain a bound on the sum of the in-band zeroes, x 2 , of S(uj' 2 ), and note that since the 
quantity S (uj' 2 )Ui(uj' 2 — uj 2 ) is an iV — 1 dimensional polynomial in uj' 2 , the sum of its zeroes 
can immediately be given as 

N-l N sr^N AN 

E2 \ ^ 2 Z^i=l \ ^ 2 ^4 /o^7\ 

Zi = - ,,2 = Z^i - — > ( 2? ) 

i=l i=l 2^1=1^1 t=l ^2 

where 

1 N 

A*n = ^E«? • (28) 

i=l 
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On denning the positive semi-definite quantity a via 

- £ = (a/) - <^ 2 ) 2 = <[^ 2 - <^ 2 >] 2 > = a , (29) 
through use of Eq. (j3J) we then obtain 

Er 2 2l 2 a 2 ^4xx(0,0) a /on\ 
Fi - ^ ] = W max - M2 = ^ max 77 • I 30 ) 

i=l ^2 M A*2 

Now in general it can be shown that the quantity A xx (0, 0)/M is related to cc^ ax according 



to 4 



UJ 



2 



2A XX (0,0)(1 + Z) 



M 

where the quantity Z is given by 



(31) 



A XX M { ] 

as summed over the even neighbors of the i = site alone. With this sum only beginning 
with the second nearest neighbors, to good approximation we can neglect the contribution 
of Z, and can thus rewrite Eq. (j3~Uj) as 

Er 2 2n ^max u /on\ 

Fi - = — — • (33) 

i=l 2 ^2 

With a being positive, we thus see that the quantity Z^l^F? — U I\ cannot be larger than 
^max/^, with the zeroes of S(u' 2 ) on average being no more than midway between adjacent 
asymptotes rather than close to the immediately adjacent higher ones. 

Further support for this result can be obtained by considering Eq. (|23j) in the infinitely 
heavy defect limit in which M' = oo, p{uj' 2 ) = — 1, corresponding to the line CD in Fig. (1). 
In this case Eq. JSJ) reduces to 

EK-^l = -^^ = -% > ( 34 ) 

with the structure of Fig. (1) then yielding 

+ {J 2 - ul) + (u'i - ul) + ... - u 2 max > E - ^] - i , (35) 



u 2 7V - 1 

max 

2 



from which the bound 



iV-l ,2 

E F* - ^] < ^ (36) 

8=1 Z 
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then follows. 

Having now established the bound on the zeroes of S(u' 2 ) which is given in Eqs. ([33)1 
and (|3T)|). we now note that for the pure mass defect case Eq. (JSJ) requires that 

yA, / 2 , , 2 1 -i- , ,' 2 , ? - Axx(0,0) A xx (0, 0) 

i=l 

where w^ ax is the largest shifted frequency. Since in the M' < M case the in-band modes 
have to lie to the left of the x\ zeroes of S(uj' 2 ), the in-band Y^=\\ jJ 'i~ {jJ T\ can then never be 
any larger than a; 2 nax /2 ~ A xx (0, 0)/M. Consequently, no matter what the value of a of Eq. 

(a quantity which varies from one crystal to the next), and no matter which particular 
crystal one might choose, if M' is less than M/2, the largest perturbed eigenfrequency ^ ax 
would than have to be larger than the pure crystal maximum u; 2 ^, and not only would it 
then lie beyond the pure crystal phonon band, once M' has been reduced enough to move 
the in-band modes as far to higher frequency as they are able to go (viz. to the zeroes of 
S(uj' 2 )), further reduction in M' would then cause w^ ax to move further and further away 
from the band maximum. Hence while the condition M' < M is a necessary condition for 
producing a mode beyond the band maximum in simple, body-centered and face-centered 
cubic crystals, it is M' < M/2 which is the sufficient one, a result previously obtained 4 by 
entirely different means 4 

IV. THE FORCE-CONSTANT CHANGE CASE 

Because of Newton's third law of motion, in the presence of force-constant changes the 
matrix V a p{l, £') of Eq. (jOJ will necessarily involve more atoms than just the one at the defect 
site and the problem is essentially intractable other than numerically if V a p(£,£') extends 
beyond the nearest neighbors of the defect. Moreover, even in the event that one restricts 
the force-constant changes to the defect and its nearest neighbors alone, for the simple cubic, 
body-centered cubic and face-centered cubic crystals the defect and its nearest neighbors 
respectively consist of a cluster of seven, nine and thirteen atoms, making V a p(£, £') 21-, 27- 
and 3 9- dimensional in those respective cases. Fortunately, because of the cubic crystal Oh 
symmetry at the defect site the V a p(£, £') matrix can be block diagonalized in the irreducible 
representations of the octahedral group, with typical decompositions 1 

r sc = Ai g + Eg + + F 2 g + 3 F iu + F 2u 

11 



Tbcc — Ai g + Eg + Fig + 2F 29 + A 2u + E u + 3Fi„ + F 2m 

T fcc = + A 2g + 2£, + 2F lg + 2F 2 , + A 2u + E u + 4F lu + 2F 2u 



(38) 



Since the displacement of the defect atom itself transforms as a 3-dimensional vector, the 
defect displacements must be located entirely in the F\ u modes, with all of the other relevant 
Oh irreducible representations being built solely out of appropriate 1 linear combinations of 
the displacements of the nearest neighbors of the defect. As such the defect displacement 
appears in 3-dimensional, 3-dimensional and 4-dimensional representations in the three re- 
spective cubic crystal cases, with the defect sector thus leading to respective 3-dimensional, 
3-dimensional and 4-dimensional blocks in Eq. (|20|) each one of which (just as in Eq. (|2ip) 
is threefold degenerate. While these blocks are still quite complicated,— in the case of central 
force-constants alone (as well as in some other specific cases such as isotropic force-constant 
crystals which we discuss in Sec. V) the determinantal condition |1 — Go^l = associated 
with the F lu block can be treated completely analytically, with both the body-centered cubic 
and face-centered cubic crystal eigenmodes being foun d 1 ^ 11 to obey precisely the same Eq. 
(|2*3*Jl as before, viz. 

1 - p{u' 2 )S{uj' 2 ) = , (39) 
save only that this time p(u/ 2 ) is given by 

A xx {0,0) 



( ax M . 2co' 2 



M' u 



2 



max 



1 



* XX M 



(40) 



For our purposes here we note that since the function S{uj' 2 ) is the same one as discussed 
previously, it still has the asymptote structure shown in Fig. (1). To determine the num- 
ber of times the function l/p(u f2 ) will intersect S(u' 2 ) requires considering each possible 
combination of mass and force- const ant changes in Eq. (J4(Jj) separately. For M' < M and 
A' xx (0,0) > A xx (0,0) the quantity l/p(^' 2 ) is everywhere positive yielding the typical curve 
AX in Fig. (2) with its iV intercepts (the point A is at +oo when M' = M), with each 
crystal eigenmode having been shifted upwards to a higher frequency. For M' > M and 
A' xx (0,0) < A xx (0, 0) the quantity l/p(u/ 2 ) is everywhere less than or equal to —1 yield- 
ing the typical curve EX in Fig. (2) with its iV intercepts (the point E is at — oo when 
M' = M), with each crystal eigenmode having been shifted downwards. For M' < M but 
A' xx (0,0) < A xx (0,0) the quantity l/p(u' 2 ) has to diverge somewhere, yielding the typical 
curve AGHX in Fig. (2) with its iV intercepts if the divergence in l/p(u/ 2 ) falls inside the 

12 




FIG. 2: Plot of S(ui' 2 ) versus to' 2 . The line AX corresponds to the value of l//?(u/ 2 ) associated 
with a defect with M' < M and A' xx (0,0) > A xx (0,0), the line EX corresponds to the l/p(u)' 2 ) 
associated with M' > M and ^4^(0,0) < A xx (0,0), the curve AGHX corresponds corresponds to 
the l/p{td' 2 ) associated with M' < M and A' xx (0, 0) < A xx (0, 0), and the curve EJKX corresponds 
corresponds to the l/p(uj' 2 ) associated with M' > M and ^4^(0,0) > A xx (0,0). 

band, 12 and analogously also yielding iV intercepts if the divergence is outside the band (not 
shown).i£ Finally, for M' > M and ^(0,0) > A xx (0,Q) there will again be intercepts 
(typical curve EJKX). 14 

As we thus see, no matter what particular values M'/M and A' xx (0,0)/A xx (0,0) might 
actually take, in all cases there are precisely intercepts, and thus precisely 3A^ eigenmodes 
in the F\ u sector. However, since this exhausts the number of degrees of freedom for the 
problem, Eq. (J2(J|) cannot generate any further eigenmode solutions. Consequently, none 
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of the determinantal conditions associated with any of the other irreducible octahedral rep- 
resentations in the cluster can yield eigenmode solutions, though just as with all the rest 
of the atoms in the crystal, the u a {€) displacements associated with these other irreducible 
representations will, via Eq. (jSJ), still participate in the Fi u mode oscillations. 

With the eigenmodes associated with the F\ u sector thus exhausting the eigenspectrum, 
it will be these modes alone which will contribute to the trace condition of Eq. (|15|). And 
with the zeroes of S(u' 2 ) still obeying Eq. for a defect with M' — M in a crystal with 
central force- const ants alone and eigenmodes which then obey 

YV 2 -^ 2 1 = yV 2 -cu 2 lW 2 -uj 2 = 2 ^ (0 ' 0) - 2 ^(°'°) ( 41 ) 

/ j 1 i % J / , L i i\ 1 max max j.j j.j ' \ I 

i=l i=l iV1 iV1 

and in-band modes which can never be shifted up beyond the immediate next zeroes of S(u' 2 ) 
(typical curve AX in Fig. (2)), we can conclude that no matter which central force- const ant 
crystal we may choose, there will definitely be a localized mode if A' xx (0, 0)/A xx (0, 0) is 
greater than 3/2, with this condition thus being sufficient to guarantee localized modes in 
nearest-neighbor central force- const ant cubic crystals when there is no change in mass.^ 

V. EXTENSION TO MORE GENERAL SITUATIONS 

To extend these results to other cases, we need to consider both non-central force- 
constants and go beyond nearest neighbors. As regards the issue of non-central force- 
constants, we note that within the nearest-neighbor approximation a complete and exact 
general relation for locating the perturbed crystal F\ u modes has recently actually been 
obtained.-^ Specifically, it was found for the body-centered cubic crystal (and thus by ex- 
tension for the face-centered cubic crystal as well since it is the reciprocal lattice of the 
body-centered cubic and we work in the harmonic approximation where momentum and 
position are treated equivalently)^ that the general determinantal condition |1 — Go^l = 
in the perturbed Fi u mode sector can be written as |1 — GqV\ = A 3 = where the nearest- 
neighbor, arbitrary force-constant A is given as 

A = £ { ttf i 1 - I 1- *] - K' 2) i 1 - w) + 4 ) • < 42) 

expressed here in terms of the two quantities p{oj' 2 ) and S(u' 2 ) which were given previously 
and two additional quantities \i and R which are needed now. To define these two additional 
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quantities one needs to introduce the pure crystal lattice Green's function combination 



R = g xy (uj'; £ = 0, £' = 222) + g xy (uj'; £ = 0, £' = 220) 

" " NM 4-' [o>?(k) - u/ 2 ] L ' 1 J 

and in terms of the interatomic potential 0(r) between two pure crystal nearest neighbors 
and the interatomic potential 0(r) between the defect and a host crystal nearest neighbor, 
define pure and impure lattice force-constants via 

/ a + (3 a a \ I a + (3 a a \ 



^(0,111) 



a 

a a + (3 a 
\ a a a + j3 / 



a a + (3 a 
\ a a a + (3 J 



4^(0,0) = -8(a + p) , A' xx (0,0) = -8(a + P) 



(44) 



where 



1 , fir) 
a = --U{r)- — 



<t>'{r) 
r 

0'(r) 



(45) 



(as evaluated at the pure crystal nearest-neighbor equilibrium separation r) , with R, fi and 
fiR then being given by 

(a + (3)((3- (3)(3a -3a + (3-@)R 



R 



fiR 



a(a — a + (3 — (3) 
2(a(3-(3a) 2 



(a + P) 2 (P - /3)(3a - 3a + /? - (3) 
2(a(3-(3a) 2 R 



(46) 



a(a + (3)(a — a + (3 — (3) 

Equation (j42j) is not only a very compact relation (it only requires knowledge of two pure 
crystal Green's functions combinations, viz. g xx {uj'\ 0, 0) and -R),— but it additionally reduces 
to the previous condition given in Eq. (pl§j) whenever [iR is zero. Now while the quantity 
\iR would vanish when (3 = (3 = 0, viz. the previously discussed central force-constant 
case, it would also vanish when a = a = 0, viz. the isotropic force-constant case, and also 
whenever a/ a and (3/(3 are equal, viz. in the mixed case in which the fractional changes 
in the central and isotropic components of the force-constants are equal to each other, with 
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A' xx (0, 0)/A xx (0, 0) = (d + (3)/(a + (3) being equal to a/a — (5/(3 in all such cases. For all of 
these cases then, the condition A' xx (0,0)/A xx (0,0) > 3/2 is sufficient to guarantee localized 
modes in nearest-neighbor crystals when there is no change in mass. 

To treat cases where the quantity \x is not zero is not as straightforward, since unlike the 
g xx (u'; 0, 0) Green's function, the R combination cannot be reduced to a sum solely over 
pure crystal eigenfrequencies as the polarization vectors cannot readily be eliminated from 
Eq. PHj) . with Eq. (|%2~j) not obviously being reducible to an expression which only involves 
the pure crystal density of states v{uj 2 ). However, despite this, it is still possible to extract a 
sufficiency condition for localized modes. Specifically, even though R does involve the pure 
crystal polarization vectors, as can be seen from Eq. (pf3*j) . its poles, and thus its asymptotes, 
are nonetheless precisely the same as those possessed by g xx (uj'; 0, 0), with the asymptotes 
of A then being none other than the asymptotes of S(uj' 2 ), i.e. none other than the ones 
exhibited in Fig. (2). To determine what happens when we take p to be non-zero then, we 
need to monitor how the intercept structure displayed in Fig. (2) gets modified as we slowly 
switch p on. 

To explicitly see what happens when we switch p on, it is instructive to first return 
to the central force-constant case, and consider the situation in which we start with some 
arbitrary A^,(0,0) which is bigger than A^O, 0) and some initial M' which is equal to 
M, and them slowly start to increase M' . We thus start with a p{uj' 2 ) which is equal to 
(2a; /2 /o;^ iax )[l — A xx (0,0)/A' xx (0,0)}, and thus with a 1/ p{uj' 2 ) which is everywhere positive 
and infinite at uj' 2 = (viz. curve AX of Fig. (3)). When we now allow M' to be just a 
little bit bigger than M, p(uj' 2 ) will now take small negative values at the smallest uj' 2 , and 
then quickly revert back to being positive again as uj' 2 is increased. Consequently, 1/ ' p(uJ 12 ) 
will now take large negative values at the smallest uj' 2 , while also quickly reverting back to 
being positive again as uj' 2 is increased (viz. curve EJKX of Fig. (3)). As can therefore 
be seen from Fig. (3), the net effect of letting M' be just a little bit bigger than M is that 
rather than increasing, the frequency of the lowest mode is instead decreased. Now since the 
trace condition of Eq. (fTB]) holds for both of these two cases, we can write a trace condition 
for the difference between the two cases, viz. 

f>f (M' > M) — uj' 2 (M' = M)\ = _ ^k|M , (47) 

i=l 

with the net effect of a small difference between M' and M on the right-hand side of Eq. 
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FIG. 3: Plot of S(ui' 2 ) versus to' 2 . The line AX corresponds to the value of l//?(u/ 2 ) associated 
with a defect with M' = M and A' xx (0,0) > ^4 XX (0,0), and the curve EJKX corresponds to the 
l/p(u' 2 ) associated with the same ^4^,(0,0) > A xx (0,0) but with M' slightly larger than M. 

(|47|) entailing a small net total shift on the left-hand side, and thus a net shift in each 
individual eigenmode of order 1/N of the shift on the right-hand side. Moreover, with the 
change in the right-hand side of Eq. ()47|) being continuous, the change on the left-hand 
side must be continuous as well. Consequently, with the transition from curve AX to curve 
EJKX needing to also be continuous, the lowest impure crystal eigenfrequency must then 
have continuously traversed the lowest lying of the pure crystal asymptotes on its way, with 
the lowest lying pure crystal eigenfrequency then being an impure crystal eigenfrequency 
at the point at which the asymptote is reached. Since S(u' 2 ) is infinite at any pure crystal 
eigenfrequency, for such a pure crystal eigenfrequency to also be a mode of the impure 
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crystal, l/p(u' 2 ) must equally be infinite at the pure crystal eigenfrequency, and it indeed 
is in this particular case since an increase in mass and an increase in force-constant act 
oppositely in piyj' 2 ), to thereby allow piyj' 2 ) to indeed vanish at the lowest pure crystal 
eigenfrequency. 

An alternate way to modify the curve AX is not to change the mass at the defect site at 
all, but to keep M' fixed at M and to instead start to allow an A' xx (0,0) which is already 
bigger than A xx (0,0) to get even bigger. Such an increase will cause all of the impure 
crystal eigenfrequencies of Fig. J3J) to increase, but will never permit any of them to ever 
reach or traverse the next immediate pure crystal asymptote since a p{u' 2 ) = (2u' 2 /u^ iax )[l — 
A XX {Q, 0)/A' xx (0, 0) which only involves force-constant changes can never vanish at any non- 
zero u/ 2 . We thus recognize two possible options as we start to vary parameters, namely that 
an impure crystal eigenmode can cross a pure crystal asymptote if p(u/ 2 ) has a zero there, 
or cannot do so if p(u/ 2 ) has no zero. If we thus picture the N asymptotes in Fig. (3) as 
dividing the uj' 2 > region into iV + 1 compartments (N — l of which lie between pure crystal 
eigenfrequencies, with the other two lying below the lowest pure crystal eigenfrequency and 
above the largest one), we see that small changes in parameters can cause the locations 
of eigenmodes to either move slightly within any given compartment or to cross into an 
adjacent one, doing so in either case without radically altering the value of the total in-band 
difference X^I 1 ^ 2 — ^f] contained on the left-hand side of Eq. (jl5j) . 

To see how this analysis pertains to non-central force- const ant crystals with non-zero p, 
we need to determine whether or not switching on p, can cause modes to change compart- 
ments. To make such a determination, we note that according to Eq. (|42|). in the event of no 
mass change at the defect site, the eigenfrequencies associated with a general force-constant 
change in a nearest-neighbor crystal are given as the solutions to 

A xx {0,0)fiR 



>2\ 



p(uj 



r2\ 



2^^(0,0) -A^(0,0)] 



4^(0,0) plR 



A xx (0,0) [1-R]_ 



(48) 



Since the pure crystal Green's function R has the same set of asymptotes as S (a/ 2 ), whenever 
S(lj' 2 ) diverges, R will do so also. However since R/(l — R) is finite at frequencies at which 
R is infinite, we see from Eq. (jlH|) that in the force-constant change case no matter what 
value we allow for p,, the perturbed eigenmodes are unable to ever leave their respective 
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frequency compartments. 

Despite this though, since R does have asymptotes, it can be expected that R would 
change sign on its way between adjacent asymptotes (though it could drop to some minimum 
and then go back to the next asymptote without ever changing sign), and thus we can 
anticipate that there will be some frequencies at which R can take an assigned value equal 
to one. At such frequencies the right-hand side of Eq. ()48j) would then diverge, and since R 
(and thus S(u' 2 )) is not itself diverging at those points, such points would not be solutions 
to Eq. (HU), but would instead be points between pure crystal eigenfrequencies at which the 
right-hand side of Eq. (|4"%|) would undergo a discontinuity. Now while such a discontinuity 
cannot take the eigenmode out of its frequency compartment, if the discontinuity is to have 
any effect at all, its only possible one would be to move the eigenmode out of the F lu mode 
sector determinant altogether. Specifically, unlike the case of a perturbation of a central 
force-constant crystal where the F\ u sector determinant accounts for all of the perturbed 
crystal eigenmodes, once we introduce non-central force-constants we have to consider the 
possibility that eigenmodes could move into the determinantal blocks of |1 — Go^l which 
are associated with the other irreducible octahedral representations of the cluster as given 
in Eq. (JHHj) . However, recalling that all such eigenmodes are also controlled by Eq. (|2U|) . 
and recalling that all the pure lattice Green's functions which appear in Eq. ()20|) have the 
common asymptote structure associated with Eq. (fTHj) and thus a compartment structure 
identical to that of S(u' 2 ) itself, we see that even if the eigenmodes do migrate into different 
irreducible Oh sectors, since the left-hand side of the trace condition of Eq. (fTHjl is a sum 
over all of the modes of the crystal and not just over those associated with the F\ u sector 
alone (a sum that thus includes the entire cluster of Eq. (JHHj) anyway), and since we are only 
making a small change on the right-hand side of Eq. (fTo^) . when the eigenmodes do move to 
other irreducible representations, they can still only move to equivalent compartments with 
the same frequency ranges as the ones they started in. As such then, independent of which 
specific sector the modes actually find themselves in after crossing any R = 1 discontinuity, 
they could only cross into the same frequency compartments they were already in. 

However, once the eigenmodes have moved out of the F lu sector, they would then be in 
some other irreducible Oh representation where the _F lu mode Eq. (}4"%j) would not apply. 
Hence in these other representations additional changes in [i (beyond the ones needed to 
first bring the force-constant dependent R of Eq. ()4(j)l to one) could perhaps then move 
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the eigenmodes into adjacent compartments. However, for that to happen the eigenmodes 
would have to cross the pure crystal asymptotes, with those particular eigenmodes then 
needing to be common eigenmodes of both the pure and impure crystals at the requisite 
value of /i. However, in order for an eigenmode to be an eigenmode at all, the eigenmode 
would need to be at a zero of the full |1 — Go^|> i.e. it would not only need to be at a 
zero of the determinantal block of some particular irreducible sector, but also to not give 
an infinite value to any of the other blocks which multiply the given determinantal block of 
interest in the full |1 — Go^l of the full crystal. However, inspection of the F\ u mode A of 
Eq. (pf2*|) shows that it is actually divergent at a pure crystal eigenmode (where it diverges as 
S(u' 2 )R), and thus once there is no asymptote crossing in the F ln mode sector, there cannot 
be any in any other irreducible representation either. Thus regardless of whether or not the 
R = 1 discontinuities take the impure crystal eigenmodes out of the F\ u sector, the modes 
cannot leave the frequency compartments they were already in when /i was zero. Hence 
if a mode already was in the compartment which lies beyond the pure crystal maximum 
eigenfrequency when the parameter [i was equal to zero, the mode will remain beyond the 
band no matter how /i might then vary. Moreover the in-band modes will remain in which 
ever compartments they had been in when /i was zero. 

We thus recognize two specific consequences of switching /z on, namely that modes have to 
stay within their respective fi = compartments, but that they are no longer prevented from 
traversing the zeroes of S(u' 2 ). Thus suppose we start off with a central force-constant crystal 
with n = and M' = M, and with ^(0,0)7^(0,0) = (a + (3) /(a + (3) = a/a = P//3 
conveniently set right at the threshold value of 3/2, and then slowly switch \i on. The in- 
band modes could now move up or down, moving to no lower than the immediately previous 
pure crystal eigenmodes (net shift of the in-band mode Y^=\\^'i ~ ^f] of down to -w^ ax ), 
or moving upwards possibly as far as the next immediate pure crystal asymptotes (net shift 
of up to a^ajj jl£ Then regardless of what value \i may actually take and regardless of which 
particular irreducible representation any specific eigenmode may actually lie in, Eq. (jjT| 
will nonetheless still hold when there is no mass change since the trace is taken over the 
entire set of irreducible Oh representations and not just over the F\ u mode. On recalling 
that = 2A XX (0, 0)/M for a nearest-neighbor crystal, we see that the necessary condition 
that wf ax -WniB be positive is given by having the in-band modes move down as far as they 
possibly can, to thereby require ^4^.(0,0) to be positive, while the sufficient condition is 
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given by having the in-band modes move up by as far as they possibly can, viz. 

w « _^ > 2^(0,0) 4^(0,0) 

"max w max — ^ M ' 

with the condition that A' xx (0,0)/A xx (0,0) be greater than two thus being sufficient to 
guarantee the generation of localized modes by force-constant changes alone in an arbitrary 
nearest-neighbor cubic crystal. 

To proceed beyond nearest neighbors is also not straightforward as no exact analog of 
Eqs. (jnni) arid is currently known. (For the case of second-nearest neighbors, the second- 
neighbor generalization of the irreducible decomposition of Eq. (JHHj) has been used to bring 
the problem to a reasonably manageable form.—) However, in the event of non-nearest 
neighbor force-constants, even though the cluster required for the matrix V a p(£', £") in Eq. 
(j2T)j) would be much larger, and even though more pure lattice Green's functions would be 
required, nonetheless every single one of these Green's functions would still have exactly the 
same asymptote and compartment structure as S(u>' 2 ). Thus if we start at the eigenmodes of 
a general nearest-neighbor crystal with non-zero fi and M' = M and slowly start to increase 
the strength of the non- nearest neighbor force- const ants (central or otherwise), we would 
again be constrained by the same set of pure crystal compartments in exactly the same way 
as before (save that there would be yet more irreducible representations of Oh to migrate to 
- with all of them also contributing to the left-hand side of the trace condition), to again 
allow us to infer that even in the most general possible case imaginable (arbitrary central 
and non-central force- const ants and arbitrary number of participating neighbors), in the 
absence of any change in mass, the condition A' xx (0, 0)/A xx (0, 0) > 2 would still be sufficient 
to guarantee localized modes.— 
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The force-constants yielded by potentials for which 0'(r) is zero when the atoms are in their 
equilibrium positions are referred to as being central, while those associated with potentials 
which obey <f>'{r)/r = <j)"(r) at equilibrium are referred to as being isotropic. 
Even though g xx (u'; 0, 0) involves a 3A-dimensional sum on both k and j, the a* J (k)cr^(k) 
term in the expression for g xx (ui'; 0,0) only permits the j sum to contribute once, since in the 
basis in which D a g(k) is diagonalized we can choose the polarization vectors <r£(k) so that only 
one polarization eigenvector has a non-zero x component. The function g xx (u'; 0, 0) thus only 
possesses N poles, poles then also possessed by g yy (io'; 0, 0) and g zz {oj'; 0, 0) as well, with the full 
3N dimensionality of the eigenfrequencies then being recovered via the existence of eigenmode 
equations for all three of u x (0, 0, 0), ity (0,0,0) and u z (0, 0, 0) in Eq. (|21|). To reflect this our 
notation here is that any eigenfrequency which appears once in an eigenmode summation from 
i = 1 to i = N will appear three times in an eigenmode summation from i = 1 to i = 3N. 
For any choice of the M'/M ratio, the curve 1/p(uj' 2 ) = M'/(M - M') is parallel to the OX 
axis in Fig. (1), lying above the line OX when M > M' and below the line CD when M' > M, 
being at ±cx> when (M - M') = 0±. 

An initial glance at Fig. (1) would appear to imply that no matter how close M' is to M, as 
long as M' is less than M there would be a solution to Eq. (|23|) beyond the band maximum, 
with the M' < M condition then being both necessary and sufficient. However if this solution 
is as close to u; max as u; max itself is to the in-band mode just below it (viz. within u; max /AQ, the 
effect of M' would then be to just slightly displace the band maximum, with the wave function 
of the mode then being a plane wave which is slightly distorted at the defect site but not at all 
localized to it; with it only being the condition that M' be less than M/2 which guarantees that 
the mode outside the band necessarily move away from the band edge (i.e. all of the decrease in 
M'/M in going from one to 1/2 could be used up entirely in moving the in-band modes to their 
allowed highest frequency maxima at the zeroes of S{oj' 2 ).) In passing we note that in order to 
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determine whether there is a bona fide mode outside the band which is localized to the defect 
site, instead of using the trace condition we can also simply solve Eq. (|23f) directly. However, in 
doing so we need to allow for the fact that beyond the band edge S(uj' 2 ) falls precipitously from 
its asymptote at to' 2 = u 2 Specifically, since S(lo' 2 ) is completely well-behaved outside the 
band and nowhere near infinite in value there, we can define its value at the band maximum as a 
continuation of its values outside the band. This is equivalent to replacing S(uj' 2 ) by its principal 
value Sp(lo' 2 ) (essentially one half of the sum of the values of S(uj' 2 ) on the two sides of the 
asymptote), to then yield an Sp^^ax) = &m which is completely finite. With S(uj' 2 ) decreasing 
monotonically beyond the band maximum, localized modes will then appear outside the band 
when M/M' > 1 + 1/S m . Since it can be shown on general grounds 4 that for any cubic crystal 
the quantity S m defined this way has to lie in the range 1 < S m < oo (with the particular value 
that S m takes in any given crystal being dependent on that crystal), the quantity 1 + 1/S m 
must lie in the range 1 < 1 + 1/SVn < 2, with M/M' > 2 thus being a sufficient condition^ for 
producing a localized mode via a mass defect in any cubic crystal. 

For the face-centered cubic crystal for instance the normalized four-dimensional basis of displace- 
ments for the AF\ U mode (ao, cti, ct2, a^) is given (in face-centered cubic crystal site notation) 
bya = u x (0,0,0), ai = (1/2 V / 2)K(1, 1, 0) + u x {l, 1, 0) + u x (l, 0, 1) + u x {l, 0, 1) + ^(1, 1, 0) + 
« a .(l,l 1 0)+« a ,(l,0,l)+ti a .(l,0,l)],a 2 = (l/2 v / 2)K(l,l,0)+n 3/ (l,l,0)+^(l,0,l)+^(l,0,l)- 
% (l,l,0)-^(l,l,0)-n,(l,0,l)-^(l,0,l)], a 3 = (1/2)^(0, 1, 1) +u x (0, 1, 1) +^(0, 1, 1) + 
u x (0, 1, 1)]. Degenerate with this basis are two others, one based on u y (0, 0, 0) and the other on 
n 2 (0,0,0). 
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On the AGHX curve the loss of an intercept outside the band is compensated for by having an 
additional intercept in the region where the l/p(u>' 2 ) divergence is located, with all modes below 
the line GH being shifted upwards, and all modes above the line GH being shifted downwards. 
After first intercepting the S{oj' 2 ) curve N times in the positive region before reaching the point 
where it diverges, on reaching this point the function 1/ p{uj' 2 ) then crosses into the negative 
region with no further intercepts then being encountered. 

On the EJKX curve the absence of any intercept at all in the region where the l/p{uj' 2 ) 
divergence is located is compensated for by having an additional intercept outside the band, 
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with all modes below the line JK being shifted downwards, and all modes above the line JK 
being shifted upwards. 

While the curve AX in Fig. (2) shows that it is necessary to have A' xx (0, 0)/A xx {0, 0) be greater 
than one in order to shift frequencies upwards, this requirement is only necessary for the produc- 
tion of localized modes by central force-constant changes alone, with it being the requirement 
that A' xx (0, 0)/A xx (0, 0) be greater than 3/2 which is then the sufficient one. (All of the increase 
in A' xx (0, 0)/A xx (0, 0) in going from one to 3/2 could be used up in moving the in-band modes 
to their allowed highest frequency maxima at the zeroes of S(ll>' 2 ).) In passing we note that as 
such, this sufficiency condition initially appears to be at variance with Eq. (|39|) . since Eq. (|39l) 
would produce a mode right at the band maximum if 2(1 — A xx (0, 0)/A' xx (0, 0)) = 1/S m and 
would thus suggest that A' xx (0, 0)/A xx (0, 0) > 2 is needed to produce a mode outside the band 
since S m can be small as one. However, while we had indicated earlier that in general S m could 
be as small as one, for pure central force-constant cubic crystals it turns out that S m is actually 
altogether larger;* making the A' xx (0,0)/A xx (0,0) > 3/2 condition the more powerful. For pure 
central force-constant crystals S m is formally found to take the value S m = oo, a value which 
however very quickly comes down to order one when there is a small non-central force-constant 
admixture - for rare gas solids for instance S m is given* by S m = 3 even though there is only a 
two percent non-central component. 

P. D. Mannheim, Arbitrary Force-Constant Changes in the Crystal Impurity Problem, 
cond-mat/0512188| (submitted to Phys. Rev. B). 

While not sought here, one can anticipate a qualitatively similar form for the determinant A in 
the simple cubic case. 

The utility of Eq. (|42|) is that while the sum over all of the F\ u mode basis neighbor sites 
would initially generate a rather large number of pure lattice Green's functions in Eq. (|2U|) . 
through repeated use of the Green's function relations of Eq. (|19|) one is, in the nearest-neighbor 
approximation, able*£ to reduce the number of needed Green's functions to a primitive set of 
just g X x(u'; 0, 0) and R. 

The presence of any analogs of the zeroes of S(uj' 2 ) in the other irreducible Oh sectors could 
prevent the eigenmodes from actually being able to go all the way to the next adjacent asymp- 
totes. 

K. Lakatos and J. A. Krumhansl, Phys. Rev. 175, 841 (1968). 
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To be slightly more precise, in relating u^ ax to 2A XX (0, 0)/M in the non-nearest neighbor case, 
one should also take into account the factor Z as given in Eq. (jHU, though we neglect this effect 
here. 
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